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ABSTRACT 

Aims. We investigate the star formation activity in a young star forming cluster embedded at the edge of the RCW 41 H n region. As 
a complementary goal, we aim to demonstrate the gain provided by wide-field adaptive optics (WFAO) instruments to study young 
clusters. 

Methods. We used deep, JHK S images from the newly commissioned Gemini-GeMS/GSAOI instrument, complemented with Spitzer 
IRAC observations, in order to study the photometric properties of the young stellar cluster. GeMS is a WFAO instrument that 
delivers almost diffraction-limited images over a field of ~2 'across. The exquisite angular resolution allows us to reach a limiting 
magnitude of / ~ 22 for 98% completeness. The combination of the IRAC photometry with our JHK S catalog is used to build color- 
color diagrams, and select young stellar object (YSO) candidates. The JHK S photometry is also used in conjunction with pre-main 
sequence evolutionary models to infer masses and ages. The K -band luminosity function is derived, and then used to build the initial 
mass function (IMF) of the cluster. 

Results. We detect the presence of 80 YSO candidates. Those YSOs are used to infer the cluster age, which is found to be in the range 
1 to 5 Myr. More precisely, we find that 1/3 of the YSOs are in a range between 3 to 5 Myr, while 2/3 of the YSO are < 3 Myr. When 
looking at the spatial distribution of these two populations, we find evidence of a potential age gradient across the field that suggests 
sequential star formation. We construct the IMF and show that we can sample the mass distribution well into the brown dwarf regime 
(down to ~ 0.01 Mq). The logarithmic mass function rises to peak at ~0.3 M©, before turning over and declining into the brown dwarf 
regime. The total cluster mass derived is estimated to be 78 + 18 Af 0 , while the ratio derived of brown dwarfs to star is 18 + 5 %. 
When comparing it with other young clusters, we find that the IMF shape of the young cluster embedded within RCW 41 is consistent 
with those of Trapezium, IC 348, or Chamaeleon I, except for the IMF peak, which happens to be at higher mass. This characteristic 
is also seen in clusters like NGC6611 or even Taurus. These results suggest that the medium-to-low mass end of the IMF possibly 
depends on environment. 

Key words. Stars: formation - circumstellar matter - ISM: bubbles - H n regions - Infrared: stars - Instrumentation: adaptive optics 
- Instrumentation: high angular resolution 


1. Introduction 


Massive stars form H n regions that expand in the surrounding 
medium. At supersonic speed, this expansion can create a layer 
of cold neutral material that is accumulated between the ion¬ 
ized and the shock fronts. This layer may become unstable and 
then collapse into fragments that form a new generation of stars 
through different physical mechanisms (see Deharveng et al. 
20101) . including the collect and collapse processes (Elmegreen 
& Lada, 197 7) or small and large scale instabilities. Observa¬ 
tions in our Galaxy suggest a possible causal link between mul¬ 
tiple generations of stars observed in star forming regions (e.g., 
Thompson et al. 2012, Preibisch et al. 120121 Dirienzo et al. 120121 
Deharveng et al. 120121 Jose et al. 120131 Samal et al. 120141) . In¬ 
terpreting these observations as evidence of triggered star for¬ 
mation requires detailed simulations, including all the possible 
sources of feedback, and several studies are currently exploring 
this aspect of star formation (e.g., Haworth et al. 120151 , Dale et 
al. l20Bl Gritschneder et al. 120141 Walch S. K. l20l4t . 


On the other hand, the observational properties of clusters 
that formed at the edge of Hu regions are poorly known (e.g., 
Krumoltz et al. 2014 , Tan et al. 120141) . mainly owing to limita¬ 
tions in angular resolution and sensitivity. For instance, the col¬ 
lapse of those massive fragments may host the formation of high- 
mass stars, as suggested by recent statistical studies, showing 
that about 30% of condensations observed at the edges of H n re¬ 
gions do host massive stars (Deharveng et al. 120101 Thompson 
et al. 120121 Simpson et al. 120121 ). On the other side of the mass 
range, young clusters contain a large number of low-mass stars, 
protostellar objects, and brown dwarfs. As predicted by mod¬ 
els, those protostellar objects and brown dwarfs will be 2 to 3 
mag brighter in young clusters than in field brown dwarfs of the 
same mass. Young clusters are then ideal targets for studying the 
nature and mass distribution of young stellar objects (YSOs), 
including brown dwarfs, and for understanding the physical pro¬ 
cess underlying their formation (e.g., Muench et al. 120031 Levine 
et al. 120061 Harayama et al. 120081 , Oliveira et al. 2 0091 Ojha et 
al. [20091 ) . 
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Young clusters usually have a heavily obscured and dense 
environment, hence deep high angular resolution in the near in¬ 
frared (NIR) is needed to resolve the stellar populations and de¬ 
tect the fainter members. Using such observational techniques, 
young cluster members can be listed and studied in detail using 
color-color (CC) diagrams and evolutionary tracks, which will 
give their age and mass distribution. The age, mass, and spa¬ 
tial distribution of stars in the cluster give important information 
about a young cluster’s evolutionary stage and put constraints on 
star formation models (e.g., Haisch et al. 2000 , Muench et al. 
2003] Levine et al. 120061 Comeron et al. 20071 Harayam a et al. 
2008] Oliveira et al. 120091 Ojha et al. 120091 DeRose et al. [2009 , 
Gennaro et al. 120111 Preibisch et al. 120111 Rochau et al. 120111 
Santos et al. 120121 Bik et al. 120141) . 

We selected a young cluster observed at the edge of the 
Galactic ionized region RCW 41 (Rodger, Campbell, Whiteoak 
1960 ) to test how wide field adaptive optics (WFAO) can help in 
reaching these goals. In particular, this cluster has been chosen 
because it has been studied in depth before, however at seeing 
limited angular resolution, in the NIR (Ortiz et al. [20071 Roman- 
Lopes et al. 120091 Santos et al. 112012k and it hosts massive young 
stars and multiple masers emission in its vicinity as signposts 
of high-mass star formation (e.g., Walsh et al. 119981 Pestalozzi 
et al. 120051 Caswell et al. 120101 Breen et al. 120111 Voronkov et 
al. 2014 ). Combining high-angular resolution, and deep NIR im¬ 
ages, we are able to study the complete populations of this young 
star cluster, well below the hydrogen-burning limit (0.08 M©). 

RCW 41 (RA = 09 / *16 m 44 5 ; DEC = -47°56'51") is a Galac¬ 
tic Hu region located in the Vela molecular ridge (VMR). Fig¬ 
ure [T] (top) shows the UKST Super-COSMOS H a emission 
(Parker et al. 2005 ). the 3.6 pm Spitzer emission, and the 5.8 
pm Spitzer emission of RCW 41. Figure [I] (bottom) shows the 
K s -bsmd emission from Soft (Son of Isaac, NIR camera mounted 
on the 3.58 meter New Technology Telescope (NTT) from La 
Silla/Chile - Moorwood et al. 1998 ). the 3.6 pm Spitzer emis¬ 
sion, and the 4.5 pm Spitzer emission. The GeMS/GSAOI field 
is shown with the contours in Figure [T] (bottom). The large scale 
morphology of the region shows strong H a emission surrounded 
by a shell of infrared emission, which is particularly visible at 
5.8 pm. The mid-infrared emission defines a ring or bubble of 
warm dust and polycyclic aromatic hydrocarbons (PAHs) that 
encompasses the ionized area. 

This structure suggests that UV radiation from massive stars 
and high-energy photoionized electrons from inside the H n re¬ 
gion heats the dust grains and excites the PAHs from the sur¬ 
rounding environment, generating the envelope type of emission. 
When looking at Fig. [T] (top), one striking element is the pres¬ 
ence of several conspicuous filaments directed almost radially 
from the cluster. This morphology has also been highlighted by 
Santos et al. ( 2012 ). where they even show that the main polar¬ 
ization direction seems to be consistent with the direction of the 
associated filament. At the edge of the H a emission, the Spitzer 
and Soft images show a dense NIR cluster, which is associated 
with the source IRAS09149 - 4743 (marked with a cross in both 
figures ). This region has been studied in detail by Ortiz et al. 
<20071) . Roman-Lopes et al. (20091 ), and Santos et al. (2012) . 
The cluster’s region is composed of two distinct substructures, 
as highlighted in Figs. [I] and [2j The main region (called main 
cluster in the following), located in the northwest part of the 
GeMS/GSAOI field, is composed of a number of objects con¬ 
centrated around a group of stars (named Objla/b in Santos et 
al. ( 2012) ) and associated with an extended emission region. The 
source Objla (Fig. |T]) has been classified as a B0 V star, and it 
is at a distance of 1.2 ± 0.12 kpc (Roman-Lopes et al. 2 009) . 
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Fig. 1 . RCW 41 region, (a) Her (blue), 3.6 pm (green), and 5.8 pm (red) 
emissions, (b) 2.12 pm (blue), 3.6 pm (green), and 4.5 pm (red). The 
green contour in fig. (a) is the field shown in fig. (b). The green contour 
in fig. (b) is the field covered by GeMS/GSAOI in the present study. 
The red cross marks the position of the NIR source IRAS09149 - 4743. 
The plus symbols show the positions of Objla and Obj2. North is up 
and east is left. 


Pirogov et al. <20071) report the detection of CS, N 2 H + , and 1.2 
mm dust continuum emissions toward IRAS09149 - 4743, with 
the dust emissions forming an almost spherical core superposed 
to the northwestern part of the stellar cluster direction. A small 
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subcluster region is located 1.1 arcmin toward the southeast of 
the main cluster’s center, and hosts Obj2 (Fig. [I]), identified as 
an 09 v star at a distance of 1.27 ± 0.13 kpc (Roman-Lopes et 
al. 200 9). which is presumably the most massive object in the 
field. Following the conclusions drawn in Santos et al. ( 20121 ). 
we adopt a cluster distance of 1.3 ± 0.2 kpc in the rest of this 
paper. 

We present images of the RCW 41 cluster corrected for at¬ 
mospheric distortions using the Gemini South Muti-Conjugate 
Ada ptive Optics System (GeMS - Rigaut et al. 1120141 Neichel et 
a!. 12014al) and obtained with the Gemini South Adaptive Optics 
Imager (GSAOI - McGregor et al. 120041 Carrasco et al. 120121) . 
GeMS is a facility instrument for the Gemini-South telescope. 
It delivers a uniform, almost diffraction-limited image quality 
at near-infrared wavelengths over an extended FoV of 2 arcmin 
across. GeMS uses five artificial laser guide stars, up to three 
natural guide stars, and multiple deformable mirrors (DMs) that 
are optically conjugated with the main turbulence layers, result¬ 
ing in a AO corrected field that is 10 to 20 times larger than 
previous AO systems (also called single-conjugate AO - SC AO). 
The GSAOI focal plane is formed by a 2 x 2 mosaic of Hawaii- 
2RG 2048 x 2048 pixel arrays with 3.0" wide gaps. Images 
are recorded in a 85" x 85" field of view with a plate scale of 
~20 mas. GeMS/GSAOI observations provide a far deeper and 
sharper view of the cluster than previously available. These new 
images allow us to reassess the stellar contents of the cluster, 
identify YSOs among its members, provide a detailed analysis 
of the mass function, and provide further insight into the low- 
mass content of RCW 41. In Sect. [2] we present the data. Section 
[3] describes the specific data reduction tools that were developed 
for this analysis. Section [4] addresses the completeness and con¬ 
tamination of the data. Section \5\ discusses the CC and color- 
magnitude (CM) diagrams. In Sect. [6] we derive the K -band lu¬ 
minosity function (KLF), and associated initial mass function 
(IMF) of the cluster. Section [7] discusses the results and Section 
[8] gives the conclusions. 


2.2. Wide field AO performance 

From the system performance point of view, what is interesting 
in Table [T]is to compare the performance during the nights of 27 
and 30. For instance, for the K s filter, the SR has been improved 
by a factor almost 3 between those two nights, while the natu¬ 
ral seeing conditions were almost identical (1.1 "for the 27, 0.95 
"for the 30). As described in Vidal et al. ( 1201 31. for both nights 
similar MCAO loops settings were used, so the difference cannot 
be explained by instrumental adjustments. In fact, the difference 
in performance can be explained by the difference in the ver¬ 
tical distribution of the turbulence over those two nights. Even 
though the integrated turbulence strength is approximately the 
same, for an MCAO system, the performance also critically de¬ 
pends on how the turbulence is distributed above the telescope 
and on how well the main turbulent layers match the optical con¬ 
jugation of the deformable mirrors (DMs). The turbulence pro¬ 
file, also called C^(h), can be evaluated from the MCAO loops 
real-time data, following the method described in Cortes et al. 
( 120121 ) . 

The normalized h) profiles for the nights of January 27 
and 30 are shown in Fig. [3] Those profiles have been normalized 
to one, and absolute profiles may be retrieved by scaling them by 
the integrated natural seeing, which in that case is almost simi¬ 
lar for the two cases. As one can see, the vertical distribution of 
the turbulence differs significantly between the two nights, since 
the profile of the 30th is more favorable to the MCAO system, 
because most of the turbulence is indeed located close to the 
DM conjugation altitude (i.e., 0 and 9km). On the other hand, 
the vertical distribution of the 27 shows a strong layer located 
at ~11km above the telescope, very likely associated with the jet 
stream (e.g., Guesalaga et al. 2014 ). and some turbulence around 
5 km, which cannot be corrected by the MCAO system. A third 
DM conjugated to 4.5km will be installed again in GeMS during 
2015 (Neichel et al. :2014a ). which should improve the perfor¬ 
mance and the stability of the MCAO correction against different 
C^(h) profiles. 


2. Presentation of the data 

2.1. GeMS/GSAOI data 

The data were obtained between January 27and February 2, 
2013 as part of program GS-2013B-SV-413 (P.I. H. Plana). This 
program is one of the 13 programs that were selected for the 
GeMS/GSAOI science verification campaign. The details of the 
observations are summarized in Table [I] Figure [2] shows the fi¬ 
nal three-color image built from the GeMS/GSAOI data. Images 
were recorded through J, H , and X s filters, with short- and long- 
exposure times to allow a wider magnitude range to be sampled 
than is possible with only a single exposure time. Each obser¬ 
vation consists of four science fields (each one dithered by 4" 
to remove gaps between the detectors) and four adjacent sky 
frames taken 2.5'away from the science field. Those sky fields 
will be used as control fields. The averaged resolution over the 
field, measured as the FWHM of the stars on single-exposure 
frames, is reported in Table [I] as is the natural seeing. The av¬ 
eraged FWHM in X s is 95 mas, while the averaged Strehl ratio 
(SR) is 14%. This is a typical performance of GeMS/GSAOI 
under suboptimal natural seeing conditions (reported as IQ85 
condition in Neichel et al. ( 12014a )). But since RCW 41 is not 
a very dense cluster, those IQ85 conditions were acceptable for 
this program. Overall, we get an improvement of the resolution 
by a factor of 4 to 6 over natural seeing and over previous NIR 
data sets (Santos et al. 2012 ). 


2.3. Spitzer-IRAC observations 

We have complemented our JHK S data with data from the In¬ 
frared Array Camera (IRAC) from the Spitzer telescope. IRAC 
has four wavelength bands centered at 3.6, 4.5, 5.8, and 8.0 pm, 
each of which having a field of view of ~5'.2 x 5'.2. The pixel 
size in all the four bands is ~ 1722. The IRAC observations of the 
region were taken on 2007 February 17 (Program ID 30734, PI: 
Donald Figer). The corrected basic calibrated data (CBCD) im¬ 
ages were downloaded from the Spitzer Space Telescope Archive 
using the Leopard package. Mosaics were built with the native 
instrument resolution of 1722 per pixel with the standard CBCDs 
using the MOPEX (Mosaicker and Point Source Extractor) soft¬ 
ware (version 18.0.1) program provided by the Spitzer Science 
Center. 

The area covered by our JHK S images is very bright at 5.8 
and 8.0 pm owing to the diffuse emission, thus we only used the 
3.6 and 4.5 pm data in the present work. Magnitudes were ex¬ 
tracted using the point response function (PRF) fitting method in 
multiframe mode and the APEX tool developed by the Spitzer 
Science Center. The standard PRF map table^] were used to fit 
variable PRFs across the image. Point sources with peak val¬ 
ues of more than 5 cr above the background were considered as 
candidate detections. Many sources detected in the nebulosity 


1 http ://s sc. Spitzer. caltech .edu/irac/calibrationfiles/psfprf/prfmap. tbl 
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Fig. 2. RCW41 false three-color GeMS/GSAOI image combining J (blue), H (green), and ^ s (red). North is up and east is left. The scale is 
indicated in the image. The location of Objla (upper right) and Obj2 (bottom left), as defined in Santos et al. (120121 ). are symbolized by orange 
stars on the image. 


appeared to be spurious. Those spurious sources were visually 
identified and deleted from the automated detection list. Fol¬ 
lowing the IRAC Data Handbook, we adopted zero points of 
280.9 and 179.7 Jy in the 3.6 and 4.5 pm bands, respectively. 
The Spitzer detections are used in Sect. |5.2| to identify the NIR- 
excess sources of the region. 


3. GeMS/GSAOI data reduction 

3.1. Data reduction 

Data was reduced using two independent methods. The first 
method used the Iraf tasks provided by Gemini, and the second 
one was based on home-made procedures developed in Yorick 
(Munro & Dubois 1 19951 ). The first steps are identical for both 
methods: (i) creation of a master flat image based on flat images 
observed each night, (ii) creation of a master sky frame based 
on the dedicated sky images, (iii) correct the science frame from 
the master flat and the master sky, as well as the detector nonlin- 
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Table 1. Observation details. FWHMs, and SRs averaged over ~ 200 stars uniformly distributed over the field and for all the individual frames. 
The corresponding standard deviations are also given. 


Date 

Filter 

Individual 
exposure time 

Number of 
frames 

<FWHM> 

(mas) 

O' FWHM 

(mas) 

<SR> 

(%) 

O ' SR 
(%) 

Natural seeing (") 
("@ 0.55/im) 

PA 

(degree) 

Jan 21 th 2013 

J 

90s 

12 

110 

20 

4 

2 

0.65 

27 


H 

80s 

12 

115 

27 

6 

3 

0.70 

27 


K s 

80s 

12 

115 

18 

10 

2 

1.10 

27 

Jan 30 th 2013 

J 

90s 

4 

150 

20 

3 

2 

1.20 

0 


J 

45s 

1 

155 

25 

2 

1.5 

1.20 

0 


K s 

80s 

14 

71 

3 

27 

3 

0.95 

0 


K s 

40s 

3 

73 

3 

27 

2 

0.90 

0 


K s 

10s 

5 

70 

2 

31 

2 

0.75 

0 

Jan 31'* 2013 

J 

90s 

8 

155 

25 

2.5 

1.5 

0.80 

0 


H 

80s 

7 

125 

20 

6 

2 

0.75 

0 

Feb 02 th 2013 

H 

80s 

8 

120 

20 

6 

2 

0.70 

0 


H 

40s 

2 

150 

15 

3.5 

1.5 

0.70 

0 


H 

10s 

10 

135 

10 

5 

2 

0.65 

0 



Altitude (km) 

Fig. 3. Vertical distribution of the turbulence (Cn2 profile), for the nights 
of January 27 (red) and 30 (black). The red profile has been shifted by 
+ 200 meters for clarity. 


earities and different gains between each detector. At that point, 
the four detector frames are reduced, but not mosaicked in one 
extension file. This last step requires applying an instrumental 
distortion correction. 

For the Iraf method, we used a high-order distortion map 
derived during GeMS/GSAOI commissioning on an HST astro- 
metric field located in the Large Magellanic Cloud (HST Pro¬ 
posal 10753, PI: Rosa Diaz-Miller, Cycle 14). The distortion 
is corrected using the program mscsetwcs in the mscred pack¬ 
age. Once the correct World Coordinate System (WCS) is prop¬ 
agated, the four detector images are combined with mscimage 
in the mscred package. The distortion map used by the Yorick 
procedures was derived from GeMS/GSAOI observations of 


NGC28£0that we correlated with HST-ACS images. Both dis¬ 
tortion maps (from Iraf and from Yorick) lead to very similar 
results, with a residual star-positioning accuracy of ~0.2". 

This is good, but not good enough for the GeMS/GSAOI de¬ 
livered image quality, which is typically around 80mas. More¬ 
over, this distortion calibration only accounts for the static in¬ 
strumental distortion. As explained in Neichel et al. ( 2014bk 
a dynamical distortion component, which depends on the NGS 
constellation and environmental factors like the telescope point¬ 
ing, is added to the final distortion budget. This dynamical dis¬ 
tortion contribution has to be corrected frame by frame, taking 
one of them as the astrometric reference. Note that ot apply¬ 
ing this second step could lead to degradinh the resolution on 
the combined image. With the Iraf data reduction package, the 
frame cross-registration is done using the program geomap im¬ 
plemented in the imcoadd package. For each image, the result¬ 
ing fit have an RMS of less than 0.3 pix in both axis, which is 
good enough to ensure no resolution degradation. The Yorick 
procedure follow the steps described in Neichel et al. ( 2014bl ). 
High-order polynomials (15 degrees of freedom per axis) are 
used to cross-register the images together. This typically leads 
to a precision of 0.1 pixel or better. 

Following the procedure described above, each individual 
image was reduced and eventually combined by filter to produce 
three long-exposure and three short-exposure reduced images. 
Both the total exposure time and the final resolution of the long- 
exposure images are summarized in Table [2] FWHM maps for 
the three bands are shown in Fig.[4]and a false three-color image 
built with the long-exposure individual image is shown in Fig. 
[2] At the cluster distance (1.3 kpc), the field covered by GSAOI 
represents -0.5 pc, with a resolution that corresponds to 130 AU. 
Looking at Fig. [4j one can clearly see a gradient in the perfor¬ 
mance across the field, the southern part getting better correction 
than the northern part of the field. Because of a technical prob¬ 
lem, this is actually due to all data taken on the night of the 27th 
having been acquired with only two NGS (the two bottom ones 
in Fig. [4]), so the north part of the field was not properly cov¬ 
ered and corrected for low-order aberrations. As a matter of fact, 
stars are more elongated owing to anisoplanatism on that region 
of the field. The effect is somehow compensated for by the data 
acquired the other nights, using all three NGS, thereby providing 


2 http://www. gemini. edu/sciop s/instruments/g s aoi/data-format- and- 
reduction?q=/node/l 1873 
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Table 2. Reduced images characteristics. FWHMs and SRs are aver- Table 3. StarFinder parameters used for the star detection and photom- 
aged over ~ 200 stars uniformly distributed over the field. Correspond- etry measurement, for each filter. The two values given for the threshold 
ing standard deviations are also given. mean that the algorithm does two iterations, with a relative threshold at 


N sigma each. 


Filter 

Total 

exposure time 

<FWHM> 

m FWHM 

<SR> 

V SR 

Filter 

Threshold 

Correlation 

J 

2160 s. 

150 

30 

2 

1 

J 

[5,2] 

0.65 

H 

2160 s. 

140 

30 

5 

1.5 

H 

[5,3] 

0.7 

K* s 

2080 s. 

95 

25 

14 

2 

K s 

[5,4] 

0.8 


better coverage of a field and better uniformity of the PSFs. The 
impact of P SFs variations on photometry accuracy is discussed 
in Sect. 13.31 

Finally, when doing a pixel-by-pixel comparison of the re¬ 
duced and combined images obtained from both reduction meth¬ 
ods (Iraf vs. Yorick), no bias is found, and the RMS of the dif¬ 
ference is lower than 5%. We therefore conclude that both data 
reduction processes are validated and that no major impact on 
the photometry analysis is to be expected from the data reduc¬ 
tion pipeline. In the following, only the Yorick images are used. 

3.2. Astrometry 

As stated above, GeMS/GSAOI images suffer from a signifi¬ 
cant field distortion, mostly owing to the instrument optical re¬ 
lay that consists of two off-axis parabolic mirrors. These dis¬ 
tortions can be compensated for when co-adding the images; 
however, the absolute WCS need to be calibrate with an external 
catalog. Astrometry with GSAOI is nontrivial because there is a 
large discrepancy between the angular resolutions of GSAOI and 
common all-sky astrometric reference catalogs. Several sources 
taken from astrometric reference catalogs are actually found to 
be multiple systems when observed with GSAOI. For that rea¬ 
son, we chose to use the Santos et al. (2 012 ) catalog, because it 
was built with the highest, previously available angular resolu¬ 
tion. We picked a set of 100 stars from the Santos et al. (12012b 
catalog, which is uniformly distributed over the GSAOI field and 
for which the GSAOI images show no close companion. The 
WCS parameters are then computed by fitting the GSAOI X/Y 
star positions, with the RA/DEC sky positions taken from the 
catalog. Only linear transformations (shift, scale, and rotation) 
are applied at that point, which leads to an averaged astrometric 
error of 80mas (maximum difference being 200mas). The av¬ 
erage precision is less than the averaged FWHM of the GSAOI 
images, so no confusion is expected. 


3.3. Flux measurement and zero-point calibration 

Flux measurement used the following two distinct methods, one 
with DAOPhot (Stetson et al. 11987 ) and the other based on 
StarFinder (Diolaiti et al. 82000 '). Before describing the steps fol¬ 
lowed in both paths, one important point to emphasize is that the 
AO PSF is fundamentally different from classical seeing-limited 
PSFs. AO PSFs usually show a complex structure, which com¬ 
bines a close-to-diffraction core, on top of rather extended wings 
(see Neichel et al. ( 2014a ) for a PSF model for GeMS). The com¬ 
plete PSFs then has a diameter that is approximately the same as 
the seeing disk, even if a large amount of the flux is concentrated 
inside an angular region of at most a few times the diffraction 
limit. For this reason, AO-photometry requires PSF-fitting algo¬ 
rithms, which provide more accurate results than regular aper¬ 
ture photometry. Moreover, the AO-PSF is not constant over 
the field or between frames. And even if MCAO instruments, 
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such as GeMS, significantly improves the uniformity of the AO- 
correction over the FoV when compared to single-conjugate AO 
(SCAO), it remains difficult to predict the PSF variations across 
the FoV. As a result, the photometric algorithms used may ac¬ 
count for PSF variation over the field, to accurately measure the 
star’s fluxes. 

DAOphot uses an analytical PSF model (e.g., Moffat), but 
allows for linear or quadratical variations across the field. Ad¬ 
ditionally, a look-up table is produced to take the deviations of 
the true PSFs from the analytical model into account. On the 
other hand, StarFinder does not use any models, but it does 
build an empirical PSF by using several stars in the image. When 
in the presence of PSF variations over the field, the choice of 
those reference stars is then one critical aspect of StarFinder 
optimization. The usual way to account for PSF variations with 
StarFinder is simply by dividing the field into subfields and 
studying each subfield individually. Of course, there is a trade¬ 
off between having enough stars per subfield, so that an accurate 
PSF model can be extracted, and to keeping the subfield as small 
as possible to keep the anisoplanatism effect to a minimum. 
The case of RCW 41 represents an almost perfect study case for 
StarFinder, because we have a large number of well-isolated 
stars uniformly distributed over the field. In that sense, RCW 41 
is not too crowded to face confusion effect (e.g., Schodel et al. 
2010 ) or too sparse for many stars to be used to build the PSF 
model. We experimented with different subfield sizes and found 
a weak dependence of the results with the number of sub-fields. 
More precisely, subfields having a size in the range between 
30"x30" and 10"xl0" were giving the same fluxes results within 
2%. For smaller fields, some subfields were affected by the small 
number of stars available to build the PSF model (<1). Farger 
fields may be affected by anisoplanatism and could suffer from 
flux errors of up to 10%. In the rest of the paper, the local PSF 
fitting is done on subfields of size 20"x20", shifted by 10" (in x 
or y) in order to introduce overlap between the frames. Each star 
will then be measured several times with different set of PSFs, 
which can be used later as a photometric error estimation. 

The other critical aspect of the StarFinder fine tuning is the 
choice of the number of iterations, the relative threshold, and 
the correlation threshold for star extraction at each iteration. We 
adjusted those parameters so that every single faint stars will be 
detected without too much contamination from the remaining 
bad or hot pixels or from structures in the background. Our three 
final reduced images (/, H , and K s ) have different characteristics 
in terms of noise and resolution, so this set of parameters has 
been optimized for each filter. Those parameters are summarized 
in Table [3] To detect as many of the faint stars as possible, we 
found that the best strategy was to set the detection threshold at 
low levels and then manually inspect each detection in order to 
remove the spurious ones. 

In the case of RCW 41, the number of stars is reasonably 
small (~ 450 / filter) so that a visual inspection and manual fil¬ 
tering of each detection is possible. Also we previously flagged 
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Fig. 4. FWHM on the reduced images. The field center is at RA = 09M6 m 44.862 5 ; DEC = -47°56'51.11". Left is J, middle is H, right is K s . The 
white stars show the NGS constellation used. Data obtained on the 27th only used the 2 bottom NGS, which explains the loss in performance 
toward the north of the field. 


all the saturated stars, which were defined assuming a minimum 
detector full well depth of 50000 e-, so that they are excluded 
from both the star extraction and photometry process. When us¬ 
ing AO, saturation may be an issue, especially for images with 
good correction. As an example, for the 80s. exposure observa¬ 
tions obtained in K s on January 30, all stars with K s < 16.5 are 
saturated in the single-frame image. To tackle this issue, short- 
exposure images were added to the observational program and 
are then used to derive the fluxes of the brightest stars. Finally, 
for four very bright stars present in the field, even the shortest ex¬ 
posure images were saturated, and we used the photometry from 
the Santos et al. (120121) catalog. The final flux catalog contains 
549 stars in K s , 483 in H, and 335 in J. Over all these detec¬ 
tions, 323 stars are commonly detected in the three filters. The 
12 stars detected in J but not in H or K s slvq either very faint 
stars, at the detection limit, or stars very close to the image bor¬ 
der where noise is higher. Combining the H and K s catalogs, 475 
are commonly detected. There are 8 stars detected in H but not 
in K s . Those stars also are on the faint end, or affected by noise 
or detector features in the ^filter. 

The next step was to convert the measured star fluxes (or in¬ 
strumental magnitude) into calibrated magnitudes. As for the as¬ 
trometry calibration, we chose to calibrate our zero-point magni¬ 
tudes based on the Santos et al. d20121 ) catalog. For that, we used 
the same set of 100 well-isolated stars. This is done both for the 
long- and short-exposure images independently. A reasonable 
match is found, with a zero-point uncertainty of ~0.1 magnitude 
for each filter. This scatter is less accurate than the zero-point un¬ 
certainty reported by Santos et al. ( 2012 ), when they calibrated 
their SofI images using 2MASS. Ss a sanity check, we also at¬ 
tempt to derive a zero point based on the 2MASS catalog. The 
mismatch in angular resolution between 2MASS and GSAOI is 
large, 2MASS resolution being ~1", while we have a factor of 
10 better, so only a tenth of well-isolated stars were identified. 
When cross-calibrating our flux measurement with 2MASS, a 
better zero-point uncertainty is found, around 0.05 magnitude 
for each filter. However, as the number of stars is reduced, sta¬ 
tistical uncertainty increases. Furthermore, those stars are on the 
bright end of the sample, where the photometric error should be 
more accurate. We also compared the fluxes measured both from 
DAOphot and StarFinder. 

Overall, very good agreement is found between the two algo¬ 
rithms, and more than 95% of the stars have instrumental mag¬ 
nitude differences between the two methods of less than 0.05 


magnitude. The remaining 5% of the stars are mainly stars ly¬ 
ing close to a bright stars, where the background and wing con¬ 
tamination from neighbors has probably not been done prop¬ 
erly. Since StarFinder was designed specifically to deal with 
AO images, we believe it is more robust to this kind of situation, 
and in the following we only use the magnitudes derived from 
StarFinder. As a conservative approach, and also to account for 
the photometric error found between the Iraf and Yorick data 
reduction processes, the zero-point accuracy is set to 0.1 mag. 


4. Photometric completeness and contamination 

4.1. Completeness and photometric uncertainties 

An important parameter is to derive the completeness limits for 
each filter. In dense stellar fields, the detectability of a source 
depends on its flux and on the local stellar density ("crowding"). 
In the case of RCW 41, detectability may also be affected by the 
diffuse background. To determine the completeness of our im¬ 
ages, we used simulated stars, embedded in the real images. This 
method somehow follows the technique described in Gennaro et 
al. (120111) . Stars were simulated based on the StarFinder empir¬ 
ical models: for each subfield location, we used the correspond¬ 
ing local PSF model. This allows to reproduce the real PSFs 
specificities and variations over the field. Those fake stars are 
introduced at random positions on each subfields. We used three 
sets of fake stars, with a total number of respectively 100, 250, 
and 500 for the whole image. We followed this strategy to make 
sure that no bias would be introduced by the number of simu¬ 
lated stars, and artificially induced crowding. Finally, for each 
filter, we simulated a range of magnitude from 13.5 to 22.0, by 
steps of 0.5 mag. Once simulated stars are added to the images, 
we ran the StarFinder detection again, using the exact same pa¬ 
rameters as the ones used for the analysis of Sect. |3.3| We can 
then estimate the number of simulated stars that are properly de¬ 
tected and compute a completeness limit. At the same time, we 
can also estimate the photometric error by computing the stan¬ 
dard deviation of the difference between the measured and sim¬ 
ulated magnitudes. 

Results for the completeness limit are presented in Fig. [5] 
(left). The values obtained, assuming 98% of completeness, are 
(. J,H,K s )i imit = (21.8, 21.5, 20.9). This is (3.2, 3.1, 3.1) magni¬ 
tudes deeper than the limiting magnitude derived by Santos et 
al. (120121) . A completeness level of 98%, such as the one chosen 
in the study, is quite conservative compared to similar studies 
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found in the literature. We used such a conservative number to 
avoid any bias towards the faint population, which is the popula¬ 
tion of main interest in this cluster. For reference, the detection 
limits are (J,H,K s ) detect = (23.0, 22.5, 21.9). 

Based on the completeness limits, one can derive the mass 
detection limit, following for instance the Baraffe et al. (120031) 
isochrones. Assuming the 1 Myr isochrones, an extinction of A^ 
= 1.0 (see Santos et al. 120121 and Sect. |4.3| ), and a distance of 
1.3 kpc (Roman-Lopes et al. ' 2009 ), we derive that the apparent 
magnitude of a M = 0.006M© star is K s = 21, which does corre¬ 
spond to our completeness limit. We can then presumably detect 
stars well below the hydrogen-burning limit of 0.08 M©, and po¬ 
tentially prove the presence of brown dwarf candidates. Pushing 
the cluster age to 5 Myr would raise the mass detection limit to 
0.015 M q , which is still well below the hydrogen-burning limit. 
We can confidently state that the completeness limit of our ob¬ 
servations probes the brown dwarf regime for any plausible age 
of the cluster. 

Results for the photometric uncertainties are presented in 
Fig. [5] (right). Those errors must be quadratically combined with 
the 0.1 mag. zero-point error derived above. This gives the fi¬ 
nal photometric accuracy. No bias was found between the 100, 
250, and 500 simulated stars cases. This confirms that crowding 
does not have a major impact on the photometric accuracy for 
RCW41. 









4.2. Magnitude histograms 

In Fig. [6j we show the star histograms for the three GSAOI fil¬ 
ters used. For reference, histograms derived from the Santos et 
al. ([201 21 catalog are overplotted, but only for the stars covering 
the same region of the cluster. For each data set, the completion 
limit derived in Sect. l4.ll is added. 

Looking at Fig. [6] the first striking element is that the 
GeMS/GSAOI data shows, for each filter, a secondary distribu¬ 
tion of star at faint magnitudes, which could not have been de¬ 
tected by the previous studies owing to their limiting magnitude. 
The nature of this second peak is discussed in Sect. |4.3| We also 
note that, even for the common region of the histogram, some 
discrepancies between the GeMS/GSAOI catalog and the Santos 
et al. (120121) exist. We looked in more detail into these differ¬ 
ences and found that all those stars were actually multiple sys¬ 
tems that were not spatially resolved by the previous SofI images. 
This is also shown by the fact that, for the common magnitude 
range, the GeMS/GSAOI histogram detects fainter stars, because 
those stars are actually resolved in fainter systems. As an exam¬ 
ple, Fig. 13 shows some representative multiple stars, which are 
taken from the H band and properly resolved by GeMS/GSAOI 
and the corresponding un-resolved SofI images. Figure [8] shows 
a close-up view of the H band central part of the main cluster, 
illustrating the gain in spatial resolution and limiting magnitude. 


4.3. Contamination 

Foreground and background stars are contaminants for studies 
of star clusters, hampering the stellar distributions of interest. 
Young cluster environments are rich with dust, potentially re¬ 
ducing the density of background stellar contamination. How¬ 
ever, in the case of deep observations close to the galactic plane, 
such as the one presented here, background contamination may 
become significant. To estimate the foreground and background 
contamination, we first used star counts from the control field, 
corrected for the surface difference between the control and the 
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Fig. 8. Close view of the center part of the main cluster around Objla 
(RA = 09 ; T6 m 43.47 5 ; DEC = -47°56 '22.86"), by GeMS/GSAOI (top 
row), and by SofI (bottom row) taken from the H band. North is up and 
east is right, and the FoV is 25"x25". The same stretch (cutting pixels 
at 1% from the minimum and 99% from the maximum) and linear scale 
has been used for the two images. 


science fields. Those star counts are shown in Fig. [6] The con¬ 
trol field observations were done without MCAO corrections, so 
the resulting resolution and completeness limit are lower than 
those of the science fields. The strategy we followed was to then 
compare the control field star counts with star counts derived 
from the Besancon galaxy model (Robin et al. 2003 ) and use 
the model star counts for the faint magnitudes. According to the 
Besancon model, foreground stars counts are negligible, and we 
only considered background stars. Background stars taken either 
from the control field or from the Besancon model must be cor¬ 
rected for the extinction in the light of sight of the cluster. We 
estimated the extinction following the same method as in Guter- 
muth et al. (12005b . by deriving the two-dimensional extinction 
maps (.Ak maps) using the //-^colors of all stars detected in both 
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Fig. 5. Left: Completeness limit estimated from simulations. Right: photometric error estimated from simulations. Red is K s , green is H , and J is 
blue. 
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Fig. 6 . Histogram of the stars for the three filters. GeMS/GSAOI data are the black histograms. Blue histograms are taken from the Santos et al. 
( 2012 ) catalog, only for number of stars per magnitude bin for the same region as the one covered by GSAOI. For each data set, the completion 
limit has been added. Dashed histogram show star counts from the Besancon model. Red histogram are the star counts from the control field. 


filters. The extinction shows a clear gradient, with higher extinc¬ 
tion values on the northwestern part of the field, the mean cluster 
extinction value being A^ uster = 1.0. For comparison, we com¬ 
puted the same extinction index for the control field. The aver¬ 
aged A k is 0.2, and it is uniformly distributed over the control 
field. The star counts derived from the Besancon model are then 
shifted toward redder magnitudes, based on the average extinc¬ 
tion difference ( 5A K ) between the cluster and the control field: 
8A K = , 4 cluster _ ^field _ 0 8 

Similar analysis was performed for J and H. The extinction 
law, for the transformation between K s and other NIR wavelength 
is done following Rieke & Lebofsky (119851) . Star counts from the 
control field, corrected for the surface difference between control 
and science field, are shown in Fig. [6} as is the star count from 
the Besancon model. For H and ^filters, both the control field 
and model agree and tend to demonstrate that star field contam¬ 
ination is fully consistent with the secondary peak of faint stars. 
For the J -band filter, a discrepancy between the field star and the 
model count appears around 18</<20, and no clear explanation 


for this difference has been found. The J -band star counts will 
then be excluded when deriving the cluster IMF (Sect. [672] ). 

A different approach to identifying a cluster member is to 
directly differentiate likely noncluster members on the basis of 
their location in a CM diagram. Santos et al. (120121) performed 
a statistical analysis of the cluster star and control field star 
distribution in the (H x [H-K s ]) and (K s x [/-^ s ]) CM dia¬ 
gram, to draw a separation criterion between cluster members 
and contaminants. A somewhat similar approach is followed in 
Harayama et al. (120081) . where they use the (/ x [/-^ s ]) CMD to 
perform a color cut to segregate cluster members alone. If we ap¬ 
ply the color-cut defined in Santos et al. 120121 we find the results 
presented in Fig. [9] Stars located on the lefthand part of the solid 
line are assumed to be field stars, not part of the cluster. Stars 
from the control field are overplotted and are all be detected as 
field stars. Field stars are mainly faint stars, which is consistent 
with what was found in Fig. [6] Finally, we also looked at the 
spatial distribution of the stars located on the secondary peak. 
As expected for field stars, those stars are uniformly distributed 
over the field, without a clear correlation with the main cluster 
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Fig. 7. Example of multiple systems taken from the H- band, resolved by GeMS/GSAOI (top row) and un-resolved by SofI (bottom row). The 
separation of the main components is, from left to right, [120, 140, 280, 1600] mas ([156, 182, 364, 2080] AU, respectively). The H -band 
magnitude measured on the first three GeMS/GSAOI images are, from left to right [13.5,14.2, 16.6], [18.4, 18.6], and [15.8, 15.9] for the multiple 
systems, while the magnitude derived from the SofI images are 13.0, 17.6, and 15.1, respectively. 



Fig. 9. (H x [. H-K s ]) CM diagram for all stars detected by 
GeMS/GSAOI. The solid line shows the statistical separation between 
cluster members (black points) and contaminating field stars (gray 
dots), as derived from Santos et al. ( 2012 ). The dotted line is an ex¬ 
trapolation of the solid line at faint magnitudes. Red triangles are stars 
from the control field. Star s flag ged with red open squares are potential 
NIR excess stars (see Sect. |5.2| ). Stars flagged with violet open squares 
are stars detected in both H and K s , but not in J. Stars flagged with blue 
open squares are also detected only in H and AT S , but they are likely to 
be field stars. 


structures. Quantitative comparison of the contamination correc¬ 


tion methods is presented in Sect. 6.1 
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5. Color diagrams and young stellar object contents 

5.1. Color-color diagram 

The CC diagram is shown in Fig.|T0| Plotted are the 259 sources 
detected commonly in JHK S , for which the photometric error 
is lower than 0.15 mag in each band. In this diagram, most of 
the sources are distributed in the reddening band of the main se¬ 
quence, which could be a mixture of reddened field stars and 
sources having warm circumstellar dust, characteristic of young 
pre-main sequence (PMS) objects (Lada & Adams 199 2). Stars 
lying below the middle reddening vector (i.e., below the MS red¬ 
dening line) are expected to be associated with NIRexcess from 
circumstellar disk. We can then use the star location on the CC 
diagram to identify the NIR-excess objects. However, since ob¬ 
servations of star forming regions are likely to be affected by 
high levels of extinction, it would make dusty objects along the 
line of sight indistinguishable from NIR excess sources. 

As an attempt to disentangle such sources, we only selected 
stars whose excess is more than lcr (where cr is the mean color 
error) from the MS reddening line and with a (J-H) color greater 
than 1.2 mag. This color cut is chosen based both on the distribu¬ 
tion of the control field sources, which mostly lie below a ( J-H) 
color of 1.2 mag (see also Fig. 5 in Santos et al. l2012 h and on the 
IRAC-identified NIR-excess sources (see Sect. |5.2| ), which are 
also found to be located above a (J-H) color of 1.2 mag. With 
these selection criterion, we may miss a few massive sources, but 
we expect to limit the degree of contamination in the final NIR- 
excess sources catalog. The 28 selected NIR-excess sources are 
marked in the CC diagram. 

5.2. I RAC identification of young stellar objects 

The circumstellar emission from the disk and envelope in the 
case of YSOs dominates at long wavelengths, where the spectral 
energy distribution (SED) significantly deviates from the pure 
photospheric emission. It is then preferable to use the longer 
wavelengths available for identifying infrared excess stars. For 
example, Haisch et al. (120001) found a fraction of NIR-excess 
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Fig. 10. ([J - H] x [H - Ks]) CC diagram for all the sources de¬ 
tected in our GeMS FOV. The solid (fire brick) and thick dashed (dark 
orange) curves represent the unreddened (bottom) and reddened (top) 
MS and giant branches (Bessell & Brett 11988) , respectively. The blue 
solid line is the locus of intrinsic T-Tauri stars (Meyer et al. 1 19971 ). The 
colors of the curves and the T-Tauri line are converted to the 2MASS 
system using the relation given in Carpenter et al. (120011 . The parallel 
dashed lines are the reddening vectors drawn from the tip of the giant 
branch (“left reddening line”), from the base of the MS branch (“mid¬ 
dle reddening line”), and form the TTS branch (“right reddening line”). 
The known O or early-B stars are shown as solid dots. The NIR-excess 
sources identified from ([J-H] x [J-[3.6]]) and ([H-K s ] x [AT s -[4.5]]) 
CC diagrams are marked as polygons (see Fig. [TlJ. The NIR-excess 
sources identified from this diagram are marked as squares. The error 
bars in the bottom right corner show average errors in the colors. 


sources of ~80% for NGC 2024 with JHKL observations, while 
a percentage of 54% is found when only JHK bands are used. 

In this work, we used IRAC observations in combination 
with our JHK S observations to reliably identify more NIR ex¬ 
cess sources than those already found in Sect. |5.1| To do so, we 
matched the 3.6 and 4.5 pm catalogs with the deep NIR cata¬ 
logs using a 172 radial matching tolerance, which corresponds to 
one pixel of the IRAC images. We then used the ([H-K s ] x [X s - 
[4.5]]) and ([J-H] x [/-[3.6]]) color diagrams to identify NIR 
excess sources (see Fig. |TT]). In these diagrams, the sources lo¬ 
cated to the right of the MS reddening vector are likely to be 
YSOs with NIR excess. Following a similar approach to the 
one discussed in Sect. [5~n NIR -excess candidates are selected 
as sources whose excess is more than lcr (where cr is the color 
error) from the MS reddening line. 

With the above approach, we identified 57 candidate YSOs with 
NIR-excess emission, five of them in common with JHK S ones 
found in the previous section. This number is probably a lower 
limit, since many sources falling in the bright nebulous regions 
were not detected in the 3.6 and 4.5 ji m bands. 

Combining the candidate YSOs identified in Sect. |5.1| (28 
sources) with the ones identified based on IRAC photometry (52 
new sources), we end up with a total of 80 NIR-excess sources 
for the area covered by our JHK S observations. 



J — [3.6] 



K-[4.5] 

Fig. 11. Top: ([J-H] x [J- [3.6]]) CC diagram of all the common sources 
detected in our GeMS FOV. The curved solid line (firebrick) is the MS 
locus. The long dashed line (red) represents the reddening vector from 
the tip of an M6 dwarf. The IR-excess sources identified from this dia¬ 
gram are marked by open polygons (see text for details). Bottom: The 
([H-K s ] x [7T s -[4.5]]) CC diagram. The curved solid line (firebrick) is 
the MS locus of late M-type dwarfs. The long dashed line (red) repre¬ 
sents the reddening vector from the tip of an M6 dwarf. The NIR-excess 
sources identified from this diagram are marked by polygons (see text 
for details). The error bars in the bottom right corner of both diagrams 
show the average errors in the colors. 


5.3. Extinction 

A first estimation of the extinction has been done in Sect. |4.3[ 
based on the method followed by Gutermuth et al. (12005 1) . In this 
section, we want to consolidate the derived extinction value by 
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estimating it directly from the CC diagram. To derive the extinc¬ 
tion of the region, we selected only the sources located within the 
MS reddening band of the ([H - K s ] x [J - H]) diagram in order 
to avoid the contribution from giant stars along the line of sight 
and the NIR-excess sources of the region. We then derived the 
extinction values of these selected sources by tracing them back 
along the reddening vector to the TTS locus, assuming that most 
of the sources are PMS in nature. We used the extinction laws 
of Rieke & Lebofsky ( [198 51), in which Ay = E(/ - i/)/0.107 and 
Ay = E (H - K)/ 0.06. The mean extinction value for the sources 
above the TTS locus was found to be Ay = 8.5 ± 2.2 mag. This 
value is consistent with the mean extinction value derived for 
the region from spectroscopic observations (Roman-Lopes et al. 
[20091 

According to the survey by Roman-Lopes et al. ( 2009 ), stars 
labeled as Objla, Obj2, and Obj4 in the ([// - K s ] x [J - H]) 
diagram belong to the most massive objects of the region with 
spectral type BOV, 09V, and B7V, respectively. The O-type and 
early B-type MS stars shows a narrow range in their intrinsic col¬ 
ors. Adopting intrinsic (/-//) colors of OB stars from Koorneef 
et al. ( 1983 ) and using the measured (J-H) photometric colors, 
we derived a visible extinction in the range of 7.2 to 8.8 mag 
with a mean ~ 8.6 mag. This is in very good agreement with the 
estimation based on the CC diagram. Th is va lue is also perfectly 
consistent with the one derived in Sect. |4.3[ of A K = 1.0 ± 0.2 
mag, considering A K /Ay = 0.112 (Rieke & Lebofsky 198 5). For 
the further analyses, we use Ay ~ 8.5 mag as the average visual 
extinction of the region. 


5.4. Age of RCW41 


The ages of young clusters are typically derived from post-main- 
sequence evolutionary tracks for the earliest members if signifi¬ 
cant evolution has occurred or by fitting the low-mass contract¬ 
ing population with theoretical PMS isochrones. Based on an 
analysis of the (H x [H-K s ]) and (K s x [/-/<f s ]) diagrams, Santos 
et al. ( 2012 ) derived an age in the range 2.5 to 5 Myr for the re¬ 
gion. Lada & Lada (12003!) note that the embedded phase of star 
cluster evolution lasts 2 to 3 Myrs, and clusters older than 5 Myrs 
are rarely associated with molecular gas. In the case of RCW 41, 
the presence of molecular gas has been confirmed (Pirogov et al. 
120071 In addition, the detection of maser sources is an indication 
of early stages of star formation in a dense circumstellar envi¬ 
ronment. These observations suggest that the cluster is young, at 
least younger than 5 Myr. 


To ascertain the age of the stellar members, we attempted a 
quantitative age determination with the help of the near-infrared 
CM diagrams. We based our age estimation on the (/ x [/-/Q]) 
CMD shown in Fig. [12] The zero-age main sequence (ZAMS) 
of Marigo et al. ( 2008 ) and PMS isochrones from the Pisa mod¬ 
els (Bell et al. I20l4b are also shown. These PMS models are 
based on an interior model of Tognelli et al. ( 2011 1 2012 ) with 
the Allard et al. (12011k Brott & Hauschildt (120051) . and Castelli 
& Kurucz ( 2010 ) bolometric corrections and additional empiri¬ 
cal corrections as discussed in Bell et al. ( 1201 41. The isochrones 
were corrected for a distance of 1.3 kpc and a visual extinction 
of 8.5 mag (see Sect. E3J.. From this CMD, it is quite appar¬ 
ent that our /-band detection limit is deep enough to detect a 
0.01 M q star at and age of 1 Myr and Ay = 8.5 mag. From Fig. 
[T2| one can notice a vertical distribution of stars, parallel to the 
reddened MS shown, but shifted toward bluer magnitudes. This 
distribution is probably due to foreground stars, indicating that 
stars belonging to RCW 41 or farther away are on its right side. 
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Fig. 12. (/ x [/-7T S ]) color-magnitude diagram of the stars in the di¬ 
rection of the RCW 41 region. Small dots show all the stars observed 
in the target field. The squares represent the NIR-excess sources identi¬ 
fied on the basis of the NIR CC diagram, whereas the polygon are the 
sources that are using NIR plus IRAC bands. The thin solid lines are 
the PMS isochrones of Pisa models (see Bell et al. [120141 for a detailed 
description) for stars of solar metallicity and ages of 1, 3, 5, and 7 Myr 
from right to left, corrected for a distance 1.3 kpc and average redden¬ 
ing A v = 8.5 mag. The reddening vectors corresponding to the 0.01, 
0.1, 0.2, 0.5, 1, 2, and 4 M 0 of the 1 Myr isochrone are drawn in dotted 
slanting lines. 


In the absence of proper motion studies or spectroscopic in¬ 
formation, the distinction of true low-mass members from the 
contaminating field stars projected along the line of sight is dif¬ 
ficult. However, since we have identified a significant number of 
likely NIR-excess candidates based on IRAC and NIR photome¬ 
try, we therefore used these NIR-excess sources to derive an ap¬ 
proximate age of the cluster. Owing to the lower sensitivity of the 
IRAC observations, and combined with the fact that low-mass 
stars do not emit significant excess emission at JHK S bands, our 
NIR-excess sample is biased toward the massive ends. But al¬ 
though the number of our NIR-excess sources is small, it should 
nonetheless reflect the age of the region. 

To determine the age of the cluster, we only used the IRAC iden¬ 
tified NIR-excess sources, because most of them are located on 
the MS reddening line of ([J-H] x [H-K s ]) diagram, suggesting 
they possess little or no circumstellar emission at these bands. 
This should limit the impact of NIR excess flux in the age de¬ 
termination. If we focus on those sources, one can see that their 
distribution in Fig.[l2]does not fall along a single isochrone, but 
instead show a scatter in age ranging from < 1 Myr to 7 Myr. 
This type of width in the evolutionary sequence of young clus¬ 
ters is common and is probably due to the combined effect of 
variable extinction, variability, and/or sources in different evolu¬ 
tionary stages. 

One interesting feature of Fig.[l2]is the potential presence of a 
group of stars falling on the 5 to 7 Myr isochrones around /~17.5 
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and with ( J-K s ) colors in the range 2.0 to 2.4 mag, which could 
be the signature of an older population. To test whether differ¬ 
ent evolutionary stages are at work in the cluster, we divided the 
NIR-excess sources into two groups: one located above the 3 
Myr isochrone and another located below the 3 Myr isochrone. 
The spatial distribution of the NIR-excess sources is shown in 
Fig. fl3] The main massive stars studied in Roman-Lopez et al. 
(120091) are also labeled. From the spatial distribution of the dif¬ 
ferent YSO populations, it appears that the "red" and probably 
younger YSOs are mainly distributed around the main cluster 
region. In Fig. [13] we overplot isophotes from the MSX band 
D emission at 14.7 gm (Egan et al. 12003b . and in the following, 
we use those contours to spatially divide the sample across the 
field. Using this spatial separation, the younger sources are also 
mainly found to be located inside the 14.7 pm emission contours 
associated with the cluster, with 85% of the "red" YSOs located 
within this contour. The "blue" and potentially older population 
of YSOs appears to be uniformly distributed over the field, ex¬ 
cept for the northern region north of the main cluster, where no 
"blue" YSOs are detected. This distributions seems to indicate 
the presence of an age gradient when going from Obj2 toward 
Obj la and beyond, where the main cluster region (around Obj la) 
would be younger than the subcluster (around Obj2) one. This 
potential age gradient has also been suggested by Santos et al. 
( 2012 ) in their analysis of the cluster evolutionary status. They 
used the PMS turn-on point as an age indicator and show that 
most (10 out of 13) of the stars that either have already reached 
the MS or are currently leaving the PMS stage are positioned 
closer to (or within) the subcluster. If this is true, then the ob¬ 
served age spread seen in the CMDs could be due to the mixed 
population of sources at the different ages present in the region. 

If we restrict the age estimate to the region located within 
the boundary of the 14.7 pm contours (i.e., within 35 arcsec ra¬ 
dius or 0.22 pc of the massive Obj la B0 star), we found that the 
majority of the sources are less than 3 Myr old. Among these, if 
we assume that all the sources falling below the 1 Myr isochrone 
have an age < 1 Myr, then the median age of the cluster would 
be < 1 Myr. It is worth noting that the PMS models are still quite 
uncertain for ages younger than about 1 Myr (see, e.g., Baraffe 
et al. 2002 ). and one should be cautious in assigning very young 
ages to PMS stars. Then, if we consider all the YSOs candidates 
of the GeMS/GSAOI FoV, our age estimate increases to ~3 Myr, 
which tends to demonstrate that the main cluster region is indeed 
younger than the rest of the field. Obj2 is presumably responsi¬ 
ble for the ionization of the evolved HII region (see Sect. |7.3| ), 
hence presumably part of the older population of the region. In 
fact, if consider only the region around Obj 2, we estimate an age 
of 4-5 Myr, but the statistical uncertainties are large, since the 
number of YSOs is low. There are only 12 YSOs among the 57 
detected with Spitzer that are lying outside of the main-cluster 
region. 

Another point that might affect the age estimate is the extinction. 
So far, we have used a constant value over the field of Ay = 8.5 
mag. The main cluster region, which is more embedded, might 
be seen through a higher extinction value, hence artificially bias¬ 
ing the age estimate. To test this possibility, we measured the ex¬ 
tinction following the method presented in Sect. |5.3| for the main 
cluster alone. Doing so, we found a mean value of Ay = 9.0 mag, 
which is very consistent with the extinction of Ay = 8.9 derived 
by Roman-Lopes et al. ( 2009 ) from their spectroscopy study. If 
we apply a higher reddening value, for instance if we use Ay = 
10 mag instead of 8.5 mag, then the age estimate within the main 
cluster region increases to < 2 Myr, which is still younger than 
the rest of the field. We therefore conclude that extinction alone 


KLF 



Fig. 14. K luminosity function built with different contamination- 
correction methods. The black histogram shows the KLF when star 
counts from the galaxy model is used, the blue dash-dot histogram 
shows the KLF when star count from the control field is used, the red 
histogram shows the KLF when the color-cut method is used to compen¬ 
sate for contamination (see Sect. |4.3| ). Vertical error bars are counting 
statistics uncertainties computed as ^fN for the black histogram. 


cannot explain the age difference seen over the field and that the 
main part of the cluster, associated with the younger YSOs, is 
probably very young, with an age < 2 Myr. 


6. Luminosity function and cluster IMF 

Stars in clusters have roughly the same age and metallicity, and 
they are located at the same distance. In addition, in the case 
of young clusters, the effects of stellar and dynamical evolution 
are minimal, then the observed present day mass function is a 
fair representation of the underlying IMF. In this section we first 
derive the cluster’s K-band luminosity function, which we use 
later to build the IMF. 


6.1. K luminosity function 


To derive the K-band luminosity function (KLF), one first needs 
to correct for field contamination and incompleteness. Incom¬ 
pleteness is not especially an issue with this data, since the lim¬ 
iting magnitudes are deep. Ho weve r, contamination is signifi¬ 
cant, as demonstrated in Sect. |4.3| and one needs to statisti¬ 
cally subtract field stars from cluster members. In Fig. [14] we 
draw the KLF when different field subtraction techniques are 
used. For each case, the correction is done for each bin by sub¬ 
tracting the contamination star count number, normalized to the 
GeMS/GSAOI field, to the global star count measured on the 
image. All histograms provide an almost similar KLF, at least 
within counting statistics uncertainties, which shows that our 
field contamination correctly removes field stars. 
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Fig. 13. Spatial location of the YSOs identified with Spitzer. Sources with red crosses are those with an age estimate. Sources with a red square 
are the YSOs with an age < 3 Myr, and blue are the YSO with an age > 3 Myr. The green diamond is a Methanol maser (Caswell et al. 12010 ). and 
the violet diamond (off the GSAOI fields) is a water maser (Breen et al. 120111) . Other masers have been detected outside the GSAOI field in the 
northwest part of the cluster. White crosses show the massive stars spectroscopically identified by Roman-Lopez et al. (120091 . Dotted lines show 
MSX isophotes at 14.7 pm (band D).The full-line circle at the edge 14.7 pm emission identify the contour of the young main cluster region. 


6.2. IMF 


To convert apparent magnitudes to absolute ones, we consider 
the distance d = 1.3 ± 0.2 kpc and a constant foreground extinc¬ 
tion of A k = 1.0 ± 0.2 mag. The mass of each star should then be 
estimated based on a mass-luminosity relation (MLR). There are 
several MLR available, covering different mass range and ages. 
Figure p~5] displays three of them, namely models from Baraffe et 
al. ( 120031 . Baraffe et al. (119981) . and Tognelli et al. (120111) . Those 
two last theoretical models have been converted into magnitudes 
using the semi-empirical bolometric corrections by Bell et al. 
(120141) . All three models are displayed with two different ages: 
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1 Myr and 2.5 Myr. As described by Muench et al. ( 2000 ). dif¬ 
ferent models are generally in good agreement, and the main un¬ 
certainties are coming from the age indetermination rather than 
model discrepancies. Based on the analysis done in Sect. |5.4| we 
have shown that an age gradient may be present in the field, with 
a younger population concentrated on the main cluster region. To 
reduce as much as possible uncertainties coming from this age 
gradient, in the following we therefore only consider the main 
cluster region for which an age of 1 ± 1 Myr has been estimated. 
The corresponding 1 Myr MLR will then be considered. 

The IMF for the RCW 41 main cluster is presented in Fig.[T6| 
There are several sources of errors that can affect the IMF shape 
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Mass-Luminosity Relation 


IMF 



Log(M) 


Fig. 15. Mass luminosity relation. Black lines indicate the Baraffe 03 
model, red the Baraffe 98, and blue Tognelli 11. Full lines denote an age 
of 1 Myr, and dashed lines denote an age of 2.5 Myr. Dash-dotted lines 
are for when uncertainties on the distance (+0.2 kpc) and on the extinc¬ 
tion (+ 0.2 mag) are summed in quadrature. The dotted lines show, for 
a given observed magnitude, the expected errors in mass. Those errors 
are used in the IMF Fig. [16] 


(e.g., Harayama et al. 200 8). and uncertainties must be treated 
carefully. The principal sources of error in this work are the 
counting statistics, which are strongly affected by the contam¬ 
ination correction at low magnitudes, hence by low-mass stars 
and by the age uncertainty, which affects the MLR relation. In 
addition, estimating the IMF from the ^f-band luminosity func¬ 
tion may be biased by the presence of NIR-excess sources. In 
Sect. |5.2| we identified 28 NIR-excess sources based on their 
JHK S colors. 

The number of NIR-excess sources located within the main 
cluster region is ten. We do not account for the IRAC excess 
sources, because they only have a negligible excess in the K s - 
band luminosity. Once corrected for contamination, the total 
number of stars is 120, so that NIR excess will convert into a 
fraction of ~8%, which is small enough so that it should not af¬ 
fect significantly the IMF. To confirm this result, in Fig. |T6} we 
plot the IMF derived when all NIR-excess sources are removed 
from the K -band sample. Both samples agree very well within 
the counting and mass uncertainties, and no systematic bias is in¬ 
troduced owing to the inclusion of NIR-excess sources. To mit¬ 
igate the effect of NIR excess sources, Ojha et al. ( 2009 ) used 
the J -band LF instead of the KLF, because it is less affected by 
circumstellar matter, compared to what is expected solely from 
the stellar photosphere. We did not try to derive a J -band LF, 
since the J -band star-counting numbers might be affected by an 
additional statistical error (see Sect. |4.3| ), but also because we 
could not apply the color-cut contaminant selection for the J- 
band. Instead, we derived the H -band luminosity function, and 
the corresponding IMF. The H -band color should be less affected 
than the K s -band by the NIR excess, and the control-field limit¬ 
ing magnitude is also deeper in H than in K s , providing a better 
constraint on the low-mass end of the IMF. The IMF derived 



Log(Mass) 

Fig. 16. IMF derived from the KLF (black solid line). The dotted line 
shows the IMF derived when NIR excess sources are removed from the 
sample. The blue dot-dashed line shows the IMF derived from the H- 
band LF. Errors in Log (Mass) are shown as horizontal lines, and they are 
estimated based on Fig. [T5] Vertical lines show the counting statistics 
error for the black histogram. The slope of the IMF’s first part, fitting 
from the black-full line histogram, gives T - -1.25 + 0.2, and is shown 
in green. The slope of the second part is T = +0.95 + 0.2, if we remove 
the peak at 0.02 M 0 from the fit, and is also shown in green. Chabrier 
et al. 12003 ) log-normal like distributions (red full line) gives a critical 
mass m c =0.35 + 0.1 M 0 and cr - 0.45 + 0.1 M©. The total mass, when 
integrating over all the bins, is 78 + 18 M©. 


from the H -band LF also agrees very well with the ^ s -band one, 
at least within the counting-statistics error bars, which confirms 
that the impact of NIR excess sources does not affect the IMF 
shape significantly in our case. 

We chose not to include a correction factor for unresolved binary. 
This choice is motivated by the fact that, with GeMS/GS AOI, we 
benefit from a high-angular resolution, which means that a sig¬ 
nificant fraction of the binary systems (those with a physical sep¬ 
aration larger than 130 AU) are actually resolved. In their anal¬ 
ysis of NGC 3603, Harayama et al. ( 200 8) found that the impact 
of unresolved binaries on the IMF slope is less than 0.04. There 
were able to resolve binaries with separations larger than 490 
AU; as a result, we can safely state that with the GeMS/GS AOI 
resolution, the impact of unresolved binaries on the IMF shape 
is negligible. Taking errors on counting statistics and errors on 
mass estimation both into account, we can derive the slopes of 
the different IMF regimes. The first regime, covering stars from 
10 to 0.1 Mq shows the typical rise in number with decreasing 
mass into subsolar regime, with Salpeter-like power-law slopes. 
The determination of the slope is sensitive to the mass and count¬ 
ing uncertainties, but also to the bin size (Maiz Apellaniz & 
Ubeda et al. 2005 ) and the mass range chosen. 

We tried different combinations of mass limits and bin se¬ 
lection (respectively 0.05, 0.1 and 0.2 M 0 ). We also tried non¬ 
regular bin sizes following the prescription by Maiz Apellaniz 
& Ubeda et al. ( 2005 ). From these experiments, we found a dif¬ 
ference in the slope of Ar = 0.2 dex. Finally, we also compared 
the slopes measured from the IMF derived from the KLF and 
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HLF, and from the NIR-excess removed sample, but all are giv¬ 
ing very similar slopes. Those results are summarized in Table]?] 
In conclusion, a conservative result would be that the first regime 
follows a power law with T = -1.25 ± 0.2. 

We also fit the subsolar IMF regime, from ~0.4 M© down to 
0.03 M q . In the subsolar regime, RCW 41 shows a steep decline 
toward fainter magnitudes and lower masses. If we fit this re¬ 
gion with a power law, we find a slope of T=0.95. If we apply a 
Chabrier et al. (t2003b log-normal like distribution to the subso¬ 
lar regime, we derive a critical mass m c =035 ±0.1 M© and cr = 
0.45 ±0.1 M q . Those values are close, but slightly higher than 
the canonical ones derived by Chabrier et al. ( 2003 ) (m c = 0.2 M© 
and cr = 0.55 Mq). When integrating over all the bins, the total 
mass is 78 ± 18 M©. All the results are summarized in Table]?] 


7. Discussion 

7.1. Comparison of the RCW41 cluster’s IMF with other 
young clusters 

We can compare the IMF derived for the young cluster embed¬ 
ded within RCW 41 with those of other well studied young clus¬ 
ters, such asIC348 and the Trapezium (Muench et al. 20031 Luh- 
man et al. 2000 ). Chamaeleon I (Luhman et al. 2007) , NGC 2024 
(Levine et al. 120061) . NGC 2264 (Sung et al. l20l0t . NGC 6611 
(Oliveira et al. 2009), NGC 7538 (Mallick et al. l20l?l) . W3 Main 
(Ojha et al. 120091 , Bik et al. 2014 ). and NGC 3603 (Harayama et 
al. 120081 . 

As for these clusters, RCW 41 has an IMF that rise in number 
with decreasing mass into the subsolar regime, with Salpeter-like 
power-law slopes. The RCW 41 slope with r=-1.25 is, however, 
less steep than for IC 348 (T = -1.53) but very consistent with 
the Trapezium one (T = -1.21) or W3 main (L= -1.1). 

Another important feature of the young cluster’s IMF is the 
turning point in mass. From the analysis of Fig. |T6] we find that 
the RCW41 IMF peaks at a mass of - 0.3 M©. Such flattening 
of the IMF around 0.5 M© has also been observed in other clus¬ 
ters (e.g., Bastian et al. 2 010. Dib S. 120141 , Luhman et al. 120071 . 
although it seems that RCW 41 is on the higher mass range than 
other similar young clusters. For instance, the Trapezium, IC 
348, Chamaeleon I and NGC 2024 all have a mode in a range 
0.1 to 0.2 M q . On the other hand, the Taurus cluster (e.g., Luh- 
mann 1120001 Briceno et al. 2 002 ) is significantly different, with 
an IMF that peaks at - 0.8 M©. The properties of RCW 41 are 
similar to those of NGC 6611, for which Oliveira et al. (120091 
derived a peak at -0.5 M©. Fitting a Chabrier log-normal distri¬ 
bution, they derived a critical mass m c =0.40 ± 0.04 M© and cr 
= 0.56 ± 0.04 Mq, which is slightly higher than but consistent 
with the values derived for RCW 41 (m c =0.35 ±0.1 M© and cr 
= 0.45 ±0.1 Mq). It is difficult to explain the differences seen 
on the characteristic mass for the different clusters. In particu¬ 
lar, we should point out that, in the mass range where the IMF 
peaks, error bars are large because of the differences in the MLR 
depending on ages. In the case of Taurus, Luhmann (120001) and 
Briceno et al. (120021) interpreted the higher mass as an unusually 
high average Jeans mass. On the other hand, Elmegreen et al. 
(120081) show that the thermal Jeans mass only weakly depends 
on environmental factors such as density, temperature, metallic- 
ity, and radiation field. Even if different clusters show significant 
differences in their IMF (e.g., Dib et al. 12014 ). it is still difficult 
to conclude on the underlying process causing those differences. 

In the substellar regime, RCW 41 shows a steep decline 
(r=0.95) toward fainter magnitudes and lower masses. Muench 
et al. (120031) derived a consistent decrease for the Trapez¬ 
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ium (r~0.7), but a much steeper decrease for IC348 (r=2.0). 
One striking element when comparing the IMFs of IC348, the 
Trapezium, and RCW 41 is the presence of a secondary peak 
around 0.02 Mq. The presence of a secondary peak in the Trapez¬ 
ium below the hydrogen-burning limit has also been observed by 
Slesnick et al. ( 2004 ). although they find it at a slightly higher 
mass (0.05 M© for Slesnick et al. vs. 0.025 Af© for Muench et 
al.). Levine et al. ( 2006 ) also report a possible secondary peak at 
a mass of 0.035 M© for NGC 2024, even though they notice that 
the significance of this peak depends on the chosen IMF bin size. 
For RCW 41, this region of the IMF is strongly affected by the 
contamination correction, and counting statistics error are large. 
Also, because our control-field limiting magnitude is lower than 
the science field, we used the calibrated galaxy model star counts 
for this mass range. Even if the control field and galaxy model 
star counts appear to be in good agreement, we cannot rule out 
any bias. For instance, Da Rio et al. ( 2012 ) analyzed the Orion 
nebulae IMF, did not find any evidence for a secondary peak, 
and interpreted this result as the potential impact of an inaccu¬ 
rate background contamination correction in the previous stud¬ 
ies. If the significance of this secondary peak is true, one inter¬ 
esting implication of such a structure in the cluster’s substellar 
IMF would be the existence of a separate formation mechanism 
for very low-mass brown dwarfs. However, as pointed out in 
Muench et al. (120031 and Lada & Lada ( 2003 ). the significance 
that should be attached to this feature depends on the accuracy 
of the adopted MLRs for subsolar objects, and it could also be 
an artifact introduced by the MLR for brown dwarfs. 

7 . 2 . Brown dwarf fraction 

An alternative way to compare different clusters’ low-mass IMF 
is to evaluate the ratio of brown dwarf to stars. Since our data are 
complete down to 0.01 M©, we should get a good estimate of the 
brown dwarf fraction for the RCW 41 cluster. We followed the 
prescription by Briceno et al. ( 2002 ). which takes the ratio of the 
number of stars for a mass range of 0.02 to 0.08 M© divided by 
the number of stars in a mass range of 0.08 to 10 M©. Doing so, 
we derived a brown dwarf fraction of 18 ± 5%. This number is 
very consistent with the fraction of 19% derived by Oliveira et al. 
(120091) forNGC 6611 and the 14% derived for IC348 by Muench 
et al. (12003k It is also reasonably consistent with the fraction of 
22% derived for the trapezium also by Muench et al. (12003b . 
It is, however, smaller than the 26% derived by Luhman et al. 
(12007b for Chamaeleon I, and lower than the fraction of 30% 
derived by Levine et al. (12006b for NGC 2024. Finally, it is sig¬ 
nificantly lower than the -43% derived for NGC 1333 by Scholz 
et al. ( 2012 ). which is actually the largest fraction observed for 
young clusters. Levine et al. ( 2006 ) suggest that the lower frac¬ 
tion seen for IC 348 could be due to the lack of photoionizing 
stars for that cluster. In fact, isolated brown dwarfs might be 
formed through the photoevaporation of accretion disks around 
prestellar cores by late O to early B main-sequence star as sug¬ 
gested by models (e.g., Whitworth[2004| Robberto et al. 12004b . 
The Trapezium, NGC 2024, and Chamaeleon I do have the pres¬ 
ence of O and/or early B stars and for those clusters, the photoe¬ 
vaporation scenario could explain the higher brown dwarfs’ frac¬ 
tion. In the case of RCW 41, the presence of late O (Obj2) and 
early B (Objla) stars has been spectroscopically confirmed. This 
is also true for NGC 6611 (Oliveira et al. 12009b . where the clus¬ 
ter hosts 03 and 05 stars. However, both clusters have asmaller 
brown dwarf fraction. For NGC 6611, Oliveira et al. ( 2009 ) sug¬ 
gest an alternative process based on gravitational fragmentation 
of infalling gas, giving rise to filamentary-like structures within 
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Table 4. Slope and parameters of the different mass functions. 


IMF 

High mass slope 
10 < M©< 0.5 

Low mass slope 
0.3 < M©< 0.03 

Log-normal fit 

1 < M©< 0.03 

Total mass 

Main cluster 

Main cluster - no NIR excess 
Main cluster - derived from HLF 

Ti = -1.25 + 0.2 
T\ = -1.21 ±0.2 

Ti = -1.18 + 0.2 

r 2 = +0.95 + 0.2 
r 2 = + 1.1 + 0.2 
r 2 = + 1.2 + 0.2 

m c =0.35 + 0.1 M©; cr = 0.45 + 0.1 Af© 
m c =0.3 + 0.1 M q \ cr = 0.4 + 0.1 M© 
m c = 0.35 + 0.1 M q \ cr - 0.4 + 0.1 M© 

78 + 18 M© 
75 + 21 M© 
77 + 20 Af© 


which lower-mass clumps would form (e.g. Bonnell et al. 120081 
Thies et al. 1 2015 ). The presence of filamentary-like structures 
(see Fig. [T}, converging toward the NIR cluster of RCW 41, 
could be an indication that such a process is also at work in 
the formation of the brown dwarfs in RCW 41. As pointed in 
Luhman et al. (12007k it is important to highlight that small sys¬ 
tematic offsets in mass estimates can result in large differences in 
the relative numbers of stars and brown dwarfs, so that looking 
for accurate differences between the different clusters might be 
difficult. In any case, we can conclude that the brown dwarf pop¬ 
ulation represents about one in the four of the RCW 41 cluster 
members, which is consistent with other nearby young clusters 
(e.g. see Kroupa et al J2013 for a review). 

7.3. RCW41 ionizing source 

The exciting star of the H n region has not been clearly identi¬ 
fied so far. According to Ortiz et al. (12007k Objla (BO v) and 
Obj2 (09 v) are the best candidates for ionization sources of 
the RCW 41 region. At 5.8 yum, the region shows a cavity at its 
center, surrounded by PAH emission in the PhotoDissociation 
Region (PDR). The H a emission is strongest in the center of 
the nebula (see Fig. [I]), implying that the ionizing source might 
be embedded within the nebula. This phenomenon has been ob¬ 
served in many bubbles associated with H n regions (e.g., Dehar- 
veng et al. [201 Oi l. Thus, we believe that the exciting star must be 
inside the 5.8 ji m cavity that surrounds the Her emission, which 
excludes the Objla source as a potential candidate. Also Objla 
is at the extreme edge of the Her emission, which makes it un¬ 
likely to be the ionizing candidate. To identify the probable ion¬ 
izing candidates of the individual Hu regions, we then used a 
JHK S catalog of the region inside the cavity and searched for 
ionizing sources in a circular area of radius ~ 1.15 pc (i.e., the 
radius of the 5.8 /rm cavity) from the center of the bubble. Using 
the extinction laws of Rieke & Lebofsky 0985 ). the observed (/- 
H ) color and M/-spectral type calibration table of Bessel & Brett 
(119881) . and a distance of 1.3 kpc, we selected luminous sources 
that have spectral types earlier than B3 MS star. We followed the 
same prescription as described in Samal et al. (120101) of reject¬ 
ing the most likely giants and foreground sources based on (/ x 
[/ - H]) and ([/ - H] x [H - /Cs 1 ]) diagrams. After these elim¬ 
inations, we are left with only one massive O-type star located 
close to the northwestern border of 5.8 gm cavity. This source 
is associated with a group of stars in its close vicinity, possi¬ 
bly part of an exciting cluster. The source corresponds to the 
"Obj2" identified by Santos et al. ( 2012 ). whose spectroscopic 
spectral type agrees with the photometric spectral derived form 
our observations. From the above discussion, we conclude that 
the Obj2 source is the most likely ionizing star of the RCW 41 
H ii region. 

7.4. Star formation toward RCW41 

Santos et al. (2012) studied the optical polarization toward the 
young cluster. Their Fig. 11 shows a color-composite image of 


RCW41 showing the optical emission from the ionized region in 
blue and the Midcourse Space experiment (MSX) band A emis¬ 
sion. MSX band A is centered at 8.2 micron and shows the emis¬ 
sion of polycyclic aromatic hydrocarbons (PAHs) molecules. 
Distribution of the R-band polarization vector is superimposed 
on their Fig. 11. For clarity we present in Fig.[l7]the same color 
composite image: the Her emission from the SuperCOSMOS H a 
survey (Parker et al. 20051 ) is shown, along with the MSX band 
A. The optical ionized region is composed of two parts: a central 
part with bright emission and a more diffuse emission located 
in the western part. The MSX band-A emission delineates the 
external bright part of the ionized region (see Fig. [17]) outlining 
the material (gas and dust) that has been collected during the ex¬ 
pansion of the ionized gas between the ionized and the shock 
fronts (see Elmegreen & Lada 11977 ). Part of this bright MSX 
8.2 micron emission is organized in filaments that converge at 
the center of the young cluster. The convergence of those fila¬ 
ments has to be proven with dynamical molecular data that are 
not available yet. However, converging flows are an important 
physical process in the formation of high-mass stars (see Kumar 
et al. 120071 Schneider et al. 201011120121) and should be consid¬ 
ered as a possible mechanism for the formation of the clump at 
the origin of this bright cluster formation and/or to the possible 
active continuous accretion of molecular material that can feed 
the cluster and allow formation of high-mass stars. Molecular- 
line observations are needed to further address what triggers the 
formation of the clump and how the filaments are active in the 
accretion process. Another question that has to be addressed is 
whether the star formation of the bright cluster has been trig¬ 
gered by the interaction with the expansion of the ionized region 
or not. A gradient of age is observed across the cluster, the bluest 
and likely ionizing sources of the region being observed on the 
lower part of the cluster, whereas young embedded sources are 
observed in the upper part, and a dense HCNO clump is observed 
farther away (see Fig. lb in Ortiz et al. 2007 ). Active high-mass 
star formation is ongoing in this dense clump (see Hill et al. 
120051 2009 ). This suggests that star formation progresses toward 
the clump and could have been triggered by the interaction of the 
ionized region with the clump. Complementary data are needed 
to confirm or disprove this scenario. 

A more diffuse optical emission extends toward the west of 
the bright central part and is delineated by fainter MSX 8.2 mi¬ 
cron emission, both filamentary and circular, that surrounds the 
external part of the ionized region. We suggest that the ionizing 
flux from the exiting stars leaks from the central part and prop¬ 
agates into the surrounding medium as observed in other bubble 
Hu regions (Deharveng et al. 2009 ). The outer part of the ion¬ 
ized region corresponds well to the photodissociation region ob¬ 
served at 8.2 micron. Leaking of the ionized emission toward the 
southwest part of the region suggests that the ionized region is 
open in this direction (lower density region) and that the ionized 
gas is coming toward the observed. We suggest that the ionized 
region is opened toward the observer and is viewed with an an¬ 
gle (beween 0 - face on and 90- edge on) on the line of sight. 
The young forming cluster appears as a bright spot of 8.2 mi- 
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Fig. 17. Color composite image of RCW41: Her emission from the 
UKST is in turquoise, MSX 8.2 micron (band A) emission in red. The 
field is 28'x 28'. North is up, east is left. The position of Objla and 
Obj2 is also shown. 


cron emission observed in the direction of the ionized gas. We 
suggest that the cluster is forming at the surface of the ionized re¬ 
gion due to the collapse of a clump that interacts with the ionized 
gas from the expanding H n region. The origin of the clump for¬ 
mation (pre-existing, converging flow, collapse from the dense 
surrounding layer of collected material) cannot be inferred from 
the existing data. 


8. Summary and conclusions 

In this paper we have presented deep, high-angular resolution 
images of a NIR cluster located at the edge of the RCW41 
Hu region. Those images have been used to explore the stel¬ 
lar content of the cluster, as well as the star formation history 
of the region. Based on these observations, our results can be 
summarized as follows: 


1. We revealed an impressive performance of the 
GeMS/GSAOI MCAO system over a large FoV. As 
part of the technical analysis, we derived that the averaged 
resolution obtained with the GeMS/GSAOI images were 150 
mas in /, 140 mas in H and 95 mas in K s , with variations of 
less than 20% over a field of almost 100 "accros. This is a 
factor of 5 to 8 times better than what can be achieved with 
seeing-limited instruments, and it shows a significant im¬ 
provement in the spatial stability of the MCAO performance 
when compared to SC AO systems. This unique performance 
demonstrates the opportunities offered by MCAO systems 
for dense regions like young clusters, in particular to 
study objects at distances unachievable before. We believe 
that efforts should be devoted to developing new analysis 
tools, specific for MCAO systems, that would reduce the 
photometry uncertainties and optimize the scientific return. 

2. Using IRAC observations to complement our deep JHK S im¬ 
ages, we were able to identify a total of 80 YSOs. Using 
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CMDs and PMS isochrones, we derived the age of these 
YSOs in the range <1 Myr to 5 Myr. More precisely, one- 
third of the YSOs exhibit an age between 3 Myr and 5 Myr, 
while two-thirds are <3 Myr. Looking at the spatial distribu¬ 
tion of these two populations, we show a potential gradient 
over the field, with the oldest YSOs mostly concentrated to¬ 
ward the subcluster region (southeast), while the youngest 
YSOs are mainly distributed around the main cluster center 
(northwest). This age gradient suggests sequential star for¬ 
mation, with the youngest part potentially being triggered by 
the first massive star generation. 

3. The exquisite resolution delivered by GeMS/GSAOI allowed 
us to characterize the stellar mass distribution down to 0.01 
M©(which is ~10 M Jup ) and to build the IMF of the stellar 
cluster. Based on the IMF, we derived a total cluster mass of 
~78 ± 18 M q . The comparison of the RCW 41 cluster IMF 
with other nearby young clusters shows that it resembles the 
Trapezium or IC 348 clusters in terms of different regimes, 
with a Salpeter-like power law at high-to-intermediate mass, 
a flattening around 0.3 M©, and a steep decrease into the sub¬ 
solar mass range. The high-to-intermediate mass regime is 
well fit by a T= -1.25 ± 0.2 slope. The IMF peak appears to 
be shifted toward high masses when compared to Trapezium 
or IC 348, and it is more consistent with NGC 661. For the 
very low-mass regime, we derived a brown-dwarf fraction of 
18 ± 5%, which is consistent with other nearby young clus¬ 
ters like Trapezium or Chamaleon I, but significantly higher 
than IC 348. Taken together, these results suggest that the 
medium-to-low mass end of the IMF depends on environ¬ 
ment. 
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